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We study the non-Markovian decoherence and disentanglement dynamics of dissipative quantum 
systems with special emphasis on non-Gaussian continuous variable systems. The dynamics are 
described by the Hu-Paz-Zhang master equation of quantum Brownian motion. The time evolution 
of the decoherence function of a single-mode superposition is compared to the concurrence of a 
two-mode entangled state. It is verified that moderate non-Markovian influences slow down the 
decay of interference fringes and quantum correlations, while strong non-Markovian effects resulting 
from an out-of-resonance bath can even accelerate the loss of coherence, compared to predictions 
of Markovian approximations. Qualitatively different scenarios including exponential, Gaussian 
or algebraic decay of the decoherence function are analyzed. It is shown that partial revivals of 
coherence can occur in case of non-Lindblad-type dynamics. 
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^ ' During the last decade quantum information and computation has been extended from discrete systems to 

— , quantum systems vifith continuous variables such as position and momentum or the amplitudes of electromagnetic 

^^ ■ field modes. This quantum information theory of continuous variable systems has received much attention in the 

\ past few years [l|, 12, l3| and has found various applications in quantum cryptography and quantum teleportation 

I— I ' [J, lal . Great advances have been made in characterizing the entanglement properties of two- mode Gaussian states 

i-Q ' by determining the necessary and sufficient criteria for their separability [6|, |7| and by developing quantitative 

entanglement measures [a Q . Non-Gaussian continuous variable states are more difficult to treat theoretically and 

to be controlled experimentally. Therefore the got much less attention in recent years. However, especially the class 

of entangled coherent states offers the possibility to apply concepts such as concurrence, first developed for discrete 

r^ ■ systems and to study the non-Markovian disentanglement dynamics of these states. 

'~^' Due to the unavoidable interaction with the environment, any pure quantum state used in some quantum 

^^ ' information process evolves into a mixed state. Thus, a realistic analysis of continuous variable quantum channels 

►^ must take decoherence and dissipation into account. Decoherence describes the environment-induced suppression 

QQ of the quantum mechanical coherence properties and interference ability. This concept is strongly related to the 

Cf^ , measurement problem 10, 11, 12, 13] and the transition from the quantum to the classical regime [14, 15, 16, 17, Ig]. 

T-H • The time scales on which these processes take place and strategies to reduce these effects are a major topic of 

C^ [ research. Thereby, the theoretical results strongly depend on the underlying dissipative dynamics and the performed 

T-H • approximations. 
T-H \ 

K^ , Within the theory of open quantum systems [13 . |20| the dissipative dynamics are mainly described by master 



^j-^ ' equations of the reduced density matrix. Initial quantum superpositions are destroyed and quantum correlations 

^ , are lost during characteristic decoherence and separability time scales. The Markovian time evolution of quantum 

k>( ' correlations of entangled two-mode continuous variable states has been examined in single-reservoir [2ll l22i] and 

rN two- reservoir models [3, [ITI, l23[ [2J|, representing noisy correlated or uncorrelated Markovian quantum channels. 

jrt • Quantum correlations are found to be better preserved in a common reservoir. Additionally the coupling to the same 

bath variables might generate new quantum correlations between the parts of the subsystem [25, 26, 27, 28, 29|. 

The underlying Born-Markov approximation assumes weak coupling between the system and the environment to 

justify a perturbative treatment and neglects short-time correlations between the system and the reservoir. This 

approach has been widely and successfully employed in the field of quantum optics [3^ where the characteristic 

time scales of the environmental correlations is much shorter compared to the internal system dynamics. Challenged 

by new experimental evidence a growing interest in non-Markovian descriptions can be observed. Very recently 

some phenomenological [3l|, [S^l and microscopic models [33, [SJ, |33, [3a| of non-Markovian quantum channels have 
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been proposed. Using the analogy between the Hilbert space of quantized electromagnetic fields and the Hilbert 
space of quantum harmonic oscillators, the Caldeira Leggett model of quantum Brownian motion [371, |38J, |39| can be 
extended to describe the entanglement dynamics of two-mode squeezed states or two-mode entangled coherent states 
as examples of Gaussian and non-Gaussian states respectively. 

Non-Markovian effects on decoherence and disentanglement dynamics become important if the decoherence time 
scale and the bath correlation time scales are comparable. This might be the case for macroscopic superposition 
since the decoherence time is reduced for increased separation in phase space. The explicit course of the decoherence 
process depends on the relation of the characteristic time scales of system r^ and environment Tb on the one hand 
and the decoherence time scale Td on the other hand. Depending on the relation of the characteristic time scales 
(including the relaxation time scale t^ four different regimes can be distinguished: 

Tb <C Ts <C T-y (Born-Mar kovian regime), (1) 

Tb ^ Ts -^ Tj (Non-Markovian regime), (2) 

Tb ^ Ts ~ T^ {Strong-coupling legiine), (3) 

Ts < Tft <^Tj (Out-of-resonance regime). (4) 

Since the decoherence time of superposed coherent states depends on the phase space separation, the initial 
preparation mainly determines the ratio of the decoherence time to the other characteristic time scales. We will focus 
on the non-Markovian regime where different decoherence scenarios exist - for both, the dynamics of single-mode 
superpositions as well as the disentanglement of two-mode entangled coherent states. 

The paper is organized as follows. In section II, we briefly review the properties of entangled coherent states 
a a special class of non-Gaussian continuous variable states and introduce relevant decoherence and entanglement 
measures. In section III, we describe the Hu-Paz-Zhang (HPZ) master equation of quantum Brownian motion which 
is the basis for studying non-Markvoian effects. We resume the extended, two-mode version of the Caldeira-Leggett 
model for single and two-reservoir models. In section IV we present and discuss the numerical results of the decoherence 
and entanglement dynamics of a single-mode superposed coherent state compared to a two-mode entangled coherent 
state. Different scenarios are analyzed, including exponential, Gaussian and algebraic decay as well as revivals of 
decoherence. It is shown that the behavior of the concurrence of an entangled two-mode state is equivalent to the 
behavior of the decoherence function of a single-mode superposition state. Finally, a brief summary is given in section 
V. 

II. NON-GAUSSIAN CONTINUOUS VARIABLE STATES 

Great advances have been made in characterizing the entanglementproperties of two-mode Gaussian states by de- 
termining the necessary and sufhcient criteria for their separability [6, 7] and by developing quantitative entanglement 
measures [Sj, jOj. Non-Gaussian continuous variable states are more difficult to treat theoretically and to be controlled 
experimentally. Therefore the got much less attention in recent years. However, especially the class of entangled 
coherent states offers the possibility to apply concepts such as concurrence, first developed for discrete systems and 
to study the non-Markovian disentanglement dynamics of these states. 

A. Entangled coherent states (ECS) 

In the following we consider entangled coherent states [40,|4l[ as a special class of non-Gaussian continuous variable 
states. An example of a multi-mode entangled coherent state is given by 

\a+)i®\a+)2®...®\a+)N+e''^\a^)i®\a^)2<Ei...®\a-)N ,^^ 

\a,a,AI) = = . (5) 

' ' V2 + 2e-2^^l"«l^cos0 ^ ^ 



which is a superposition of coherent states \a±)i — I){±ao)\0)i with displacement parameter laoP = ^pV -I- ^jS^ and 

ctq as width of a minimum uncertainty wave packet displaced at qo, po (mean values) in phase space. In particular, 
we will focus on the two-mode version 



which is also known as quasi-Bell-state [3J, |42, |43| due to its similarity to the well-known Bell states as maximally 
entangled two-qubit states. In addition to a violation of the Bell-inequalities '4411 these states show further nonclassical 
properties, such as sub-Poissonian statistics, squeezing and a violation of the Cauchy-Schwarz inequaltiy 45] . Con- 



trary to entan gled single Fock-states [46j ECS are superpositions of non-orthogonal states [47[ with some remarkable 

properties [11,51, [13]. 

If just a single mode is considered, we get the well known cat-state superposition of two coherent states 

, \a+)+e''>\a-) 

'^"^ " V2 + 2e-2|"oPcos0 ^' 

(which is of course not entangled). Its density matrix can be separated into two parts, 

_ \a+){a+\ + \a-){a-\ \a-){a+\ + \a+){a-\ , . 

^^ V ' " V ' 

PclW P/(0) 

a classical part pd and a part pi describing the quantum mechanical ability to interfere (here = Q). The corresponding 
phase space representation in form of the Wigner function 

e-^V^-g-^^^-g/^- (^,.g/2.g,,,gg2P+eoshM^ (9) 

2 

(for Jao]^ — -pr) has two characteristic peaks at g = ±go and an oscillating pattern in between with partially negative 
values of the Wigner function. It is one of the most common examples for studying decoherence for continuous 
variable states. In this paper we want to compare the decoherence properties of these state within a non-Markovian 
description to the disentanglement process of an two-mode entangled coherent state. 

B. Decoherence and Disentanglement measures 

1. Decoherence of a single-mode superposition state 

There are different measures of decoherence in phase space, in position and momentum space or Fock space represen- 
tation. Decoherence time scales based on different measures can deviate from each other. In our view, phase space 
measures give the better picture of the underlying non-Markovian process compared to a treatment in position space 
alone (i.a. the so-called attenuation factor or the fringe visibility function). Decoherence in phase space can be de- 
scribed by the time evolution of the purity p{t) = Tr[p2(i)] or the so called quantumness S(i) = J dqdp\W{q,p, t)\ — l 
which measures the phase space volume of the negative part of the Wigner function. While the first mixes decoherence 
and relaxation phenomena and the latter is difficult to calculate, the norm /x/(t) = Tr[pj(t)] of the interference part 
seems to be an adequate measure for decoherence [51|. It can be calculated from the oscillating part Wi{q,p^ t) of the 
Wigner function by 



pi{t) = {2nhY dq dpWfiq,p,t). (10) 

J —CO J —CO 

The non-Markovian time evolution of the Wigner function is governed by a quantum Fokker-Planck equation given 
in the third section. 

2. Disentanglement of a two-mode entangled coherent state 

To study the entanglement properties of a non-Gaussian continuous variable state we use the concept of con- 
currence j52l |53|. For the density matrix pi2 of a pair of qubits the concurrence is defined as C12 = 
max{Ai — A2 — A3 — A4,0}, where Ai > A2 > A3 > A4 are the square roots of the eigenvalues of the spin-flipped 
state Q12 = P12 (cTy (g) ay) /9*2 i^y "^ '^y) with Pauli matrix ay. The concurrence varies from C12 = for disentangled 
qubits 1 and 2 to C12 — 1 for maximally entangled states. To determine the pairwise entanglement of a A^-mode 



system in analogy to the discrete A^-qubit system the reduced density matrix pi2 — Tr2.4,..Ar {|a, 0, A^)(a, 0, A^|} of 
two modes 1 and 2 has to be considered (all reduced density matrices pki of any two modes k and I are identical). 
The trick is to represent the density matrix 



1 



Pl2 



{|a+)|a_)(a+|(a_| + |a_)|a+)(a_|(a+| + 






(11) 



in an orthogonal basis {|0), |1)} with |0) = |a+) and |1) = — 



l"+> 



( 2p^{l+qcoiie) pM{l + qe'^) pAf(l + ge-*^) 0\ 

pM{l + qe-'^) IVP M^qe-''^ 

pAf(l + ge*^) APqe''^ M^ 

oy 



(12) 



P12 - K 

\ 

with parameters p = e^^'"' , A^'o = (2 + 2p^cos0)^^, q = p^^^ and M = -^/l — p^ [4j|, By determining the four 
eigenvalues Ai = NqAP{1 + q), X2 = NqM'^{1 — q), and A3 = A4 = the concurrence of the reduced state pi2 can be 
calculated to 



Cl2 = 



APq 



g-2JV|ap(l_g4|ap) 



1 + p^ COS 



1 



,-2JV|qP 



cos6' 



(13) 



In case of a two- mode state {N = 2) we have pure state entanglement. The entanglement dynamics are then described 
by the time evolution of the parameters p{t), q(t) and M{t). A specific model will be introduced in the next section. 

III. NON-MARKOVIAN DYNAMICS OF CONTINUOUS VARIABLE STATES 
A. HPZ master equation of quantum Brownian motion 



Non-Markovian effects are discussed here on the basis of the Caldeira-Leggett model of quantum Brownian motion 
[371 [sil . [53 often referred to as independent-oscillator- model [3a, \30\ ■ It is a system plus reservoir model where the 
total Hamiltonian H = Hg + Hi, + Hint consists of three parts, with Hg as Hamiltonian of the subsystem which 
interacts via the Hamiltonian iJint with a bath that is described by a collection of a large number of harmonic 
oscillators H^ = J2i fi^i{b^b -1-1). In detail the Hamiltonian of the Caldeira Leggett model is given by 



H 



L 

2m 



N 



l^+^w+i: 



2mi 



Xi 



Ciq 



(14) 



where q and p are the Heisenberg-operators of coordinate and momenta of the Brownian oscillator moving in an 
harmonic potential V{q) = -^mbj^q^ and coupled to a bath of N independent harmonic oscillators with variables Xi, 
Pi and frequencies LOi. The bath is characterized by its spectral density 



J{uj) 



N 



^-^ 2muj. 



■(5(W - LUi) 



7wr2 



r2- 



(15) 



The interaction is bilinear in the coordinates q and Xi of the subsystem and the bath respectively with coupling 
parameters c^. The self- interaction term (proportional to q^) in the Hamiltonian 



Hu 



E 



-CiX^q - 



2miUjf 



(16) 



renormalizes the oscillator potential to ensure that the observable frequency is close to bare oscillator frequency coq. 
From influence functional path integral techniques Hu, Paz, Zhang have derived the master equation 



p = - [Hs,p\ + ^./-' [q^p] + ^[9, {p, p}] 



2ih 



2ih 



^hhp]] 



—j^hhp]], 



(17) 



with [ . ] and { . } denoting commutator and anti-commutator respectively. This master equation is vahd for arbitrary 
couphng and temperature. The non-Markovian character is contained in the time-dependent coefficients which read 
in expansion up to the second order in the system-bath couphng constant }39j : 

D„p{t) = -^ f dt'K{t') sin LJot' *^' mjJoo)(q^)-^, (20) 

mujo Jq m 

Dp{t) = / dt' K{t') cos Luot' *^ (p^)7p(oo), (2f) 

Jq 

where L{t) — i{[r]{t),r]{0)]) and K{t) = 5({'7(0,?7(0)}) are connected to the spectral density (fTS]) by 

L{t) = - / dujj{u)s\\\u}t, (22) 

TT Jo 

ft /"^ f 

K{t) = - / duo J {u)coih{- Phw) cos ujt. (23) 

Ti" Jo 2 

ii'(i) is the correlation function of the quantum noise term -q resulting from averaging over the initial thermal bath 
distribution. The exact expressions of the HPZ-coefhcients are related to the Green's functions of the correspond- 
ing quantum Langevin equations |57l . |58| . The entanglement properties of the joint state of the oscillator and its 
environment have been studied in ref. [59|. 

B. Characteristic time scales 

The characteristics of the decoherence process depend on the relation between the different time scales of the 
systems on the one hand and the decoherence time scale on the other hand. The characteristic time scales of the 
system - the internal system dynamics Tg, the relaxation time scale r^ and the bath correlation time scale T(, - are 
determined by the coefficients of the HPZ master equation (fT7|) and can be approximated by 



(c.2 + <5r)2)-i/2 = 






^o\ (24) 



-. «7- = 7(^l + ^j - 7- (25) 

n ~ min{r-\ (3h} - r-^ (26) 

The decoherence time scale Td for entangled coherent states is mainly governed by the phase space separation in form 
of the parameter \a\. 

C. Secular approximation of the HPZ master equation 

Performing a secular approximation of the HPZ master equation (fT7|) by averaging over the rapidly oscillating terms 
of the time-dependent coefficients (IT51) - (PT|) (which is equivalent to a rotating wave approximation after tracing over 
the environment without neglecting the counter-rotating terms) one gets the following approximated master equation 

[13, [Ml m, [el: 



p — —iujo[a^a, p] H — [2apa' — a'ap — pa^a\ -\ — — [2a' pa — aa' p — paa''\ . (27) 

The form of this equation is similar to the quantum optical master equation of the damped harmonic oscillator in 
Lindblad form, with the only difference that the coefficients 7|.| appearing in the master equation are time-dependent. 



The connection to the HPZ-coefHcients (TTSl) and (l?T|) is given by 

7X W ^(Ml + m ^>^' ^^ (. + 1) , (28) 

'-^^ ^ \hmuJo 2 J muj^ + T^^ ' ^ ' 

In the limit i ^ F"^ the reach the corresponding Markovian values adjusted by a factor r^/(r^ + Wg) ~ 1- As long 
as the coefficients 7^ are positive for all times the equation (P7)) is of Lindblad type. However, not every master 
equation of Lindblad type form does necessarily fulfill the semigroup property [6l| . For certain parameter ranges the 
coefficients can become negative and the dynamical evolution is of non-Lindblad type [60l | . 

D. Two-reservoir model 

The dynamics of two identical, not directly interacting modes (with coordinates and momenta qj,Pj, j = 1,2) in 
two uncorrelated reservoirs is modeled by the interaction Hamiltonian 



oc 



^int = -<?1 ^ CiXi - 52 ^ CiX'r (30) 



i=l 



with (a;^x° + x^x^) = Vi, j. The master equation of the reduced density matrix is then given by the sum of the 
master equations of two single modes [3j] : 






2ihm 2ih 



^fe,fe,p}] + ^fe,b,P]]-^fe,fc,p]]|. (31) 



The time dependent coefficients are given by 7^ (i) = ujq + Sil^ (t) — 'yT/m and eq. (fT5|) to (I^T|) . The time evolution 
of a (Gaussian) two-mode squeezed state in two uncorrelated non-Markovian channels has been studied recently in 
ref. [231 (while we focus on the evolution of non-Gaussian states). The authors derived a non-Markovian separability 
function which shows oscillations in case of an artificial out of resonance bath with T <ti luq. In this two- reservoir model 
the initial entanglement is always completely lost and both modes are finally uncorrelated (even at zero temperature 
while Ti — > cx)). If the secular approximation is applied, the non-Markovian dynamics can be described by the master 
equation 



dt 



- -L— Y^ ^2ajpa] - a]ajp - pa]aj + -L— ^ \2a]paj - aja]p - pa^a] , (32) 



with time dependent coefficients given in eq. (pS)) . It should be noticed that for n — > the coefficient 71 (i) is different 
from zero for short times and vanishes for times t 3> F^^. 

E. Solution of the model 

In this section we present the solution of the HPZ master equation for the relevant initial states. 

1. Decoherence function 

Given the characteristic function xiiVj ^i 0) of an initial preparation, the time-dependent Wigner function Wi{q,p, t) 
can be calculated from [64||: 

1 /"OO /"OO 

Wi{q,p,t) = I dr] dvxi{'r]t,i't,0) exp[i{r]p + i^q)/h], (33) 

{2TThy j_^ j^^ 



where xi is identical to the initial characteristic function by substituting r/t — f{t)ri + —f{t)v, vt = 'n^I{t)v + f{t)v 
and multiplying a Gaussian factor, 



XiiVu J^t, 0) = X/ f{t)ri + -.f{t)v, mf{t)r^ + f{t)v- exp 



^ — {Kpit)rf + 2K,p{t)T^iy + K.ity) 



(34) 



The coefficients Kp^q{t) and Kqp{t) as well as the correlations {q^{t)), {p'^{t)) and {{q,p}{t)) are obtained as solution 



of the corresponding quantum Langevin equation with Greens function f(t). With a — ^%-/(i), b 
-^0 = (1 + e^o''^o)~^ the norm Tr[pj] of the interference part finally reads 



M/(0 = Mo(0^o + Mo(i)^o exp 



Mo(09o 



h?7 






r/(i) and 



(35) 



(normalized to /^/(O) = 1). This is our decoherence function. 



2. Concurrence 



The concurrence of a two- mode entangled coherent state under Markovian evolution in a zero temperature environment 
with 7x = 7 and 7| = can be calculated by introducing a time dependent orthogonal basis. 



|0(i))^|a+(i)), \l{t))^ 






(36) 



,-2e-"|Qp 



with \oi±{t)) = I ±aoe ^*/^). This leads to a time dependent density matrix pi2{t) with coefficients p{t) = 

q{t) — e^^'^^^'^ ^ ^l"l , and M{t) = \J\ — p'^{t). Under non-Markovian evolution the constant coefficients 7|, 7| have 

to be replaced by the time dependent coefficients ^i [t] , 7| (i) which leads to the substitution 



7t^ rp(i) = / dsjp(s) and A(i) -^ AM) = — / dse^'^'^DJs). 

Jo 'I'muJo Jo 

Thus, we receive the time dependent coefficients of the density matrix 
p{t) = e-^^"^"'"!"!', q{t) = e-4Ap(t)l"P^ M{t) - ^l~p^{t), 
From the corresponding eigenvalues Xi{t), i — 1, ...,A the concurrence can be calculated 



(37) 



(38) 



CMt) = N'M'{t)q{t) = 



,-4Ap(t)|a|^ 



1 _e-4|a|2(,Qg 

and is compared to the time-dependent decoherence function of a single-mode superposition state in section IV. C 



(39) 



IV. DECOHERENCE AND DISENTANGLEMENT SCENARIOS IN THE NON-MARKOVIAN REGIME 

A. Exponential and Gaussian decay 



For given system and bath parameters the decoherence time scale varies with the initial phase space separation 
|q;|. The parameter go = ^l? ^ ?'! therefore allows to determine the relation between the system dynamics t^ and the 
decoherence timescale r^. In the limits t^ 3> t^ and r^ <?C t^ it is possible to derive the decoherence time analytically 
by considering just the leading terms in the HPZ master equation p7|) . 
For Tg <^ Td <S^ T-y the diffusion coefficient Dp[t) dominates the dynamics thus having 



Mr^>r,(0 



exp 



'2^Dp{oo)t 



n \ 2 



exp 



4gi 



dt'K{t') cos ujot' I t 



(40) 



where the stationary value Dp{oo) = jp{oo){p'^{oo)) with (p^(oo)) is received from the solution of the corresponding 
Langevin equation for times t « r~^ <C Td- For 7 <C mujo and T ^ ujq we have Dp{oo) = ^ ,"° coth(2/3ftti;o), thus for 
/cT 3> ftwo receiving the decoherence function 



t^Ta>rAt) « exp 



47^0 go 
2n 



coth(-/3fixi;o)i 



kT:S>huo 



exp 



47fcTgg 
ft2 



(41) 



Decoherence is dominated by thermal fluctuations and the decoherence function decreases exponentially with charac- 
teristic time scale 

ft2 



Mrd>T,(0 "" exp[-t/rrf] with r^ 



ijkT 



% 



(42) 



and scales inverse quadratically with the displacement go [65( ■ 

In the limit Td -^ Tg ^ r^ decoherence takes place on time scales t <C Wq and is dominated by vacuum fluctuations 
and system-bath interaction. The decoherence factor can be derived from influence functional path integrals methods 
and is given by 

2 



Mr^<r,(t) 



exp 



49| 

ft2 



dt' / dt"K{t' - t") 
Jo 



exp 






1 1 

dioJ{uj) coi\i{— pfvjj) — 



cos(a;t) 



(43) 



The decay of the decoherence function becomes more and more non-exponential when the parameter q^ is increased 
and finally becomes Gaussian for Td -^ tij <^ Ts <^ T-y with 



A*Trf<r,(i) ~ cxp[-(t/Td) ] with 



Td = 



2^/im 



_i kT-»nwo 

% — ' 



2y/jrkT 



%\ 



(44) 



where the dissipation-fluctuation theorem K{t) = kTj{t) = kT^Te"^*^ has been applied. In this case the decoherence 
time scales just linearly with the inverse initial separation q^ and depends explicitly on the bath cut off frequency 
r. Determining the decoherence time by just considering the diffusion term leads to an underestimation of the 
decoherence rate. Thus, the decoherence scenarios for mesoscopic (ts <C Td <C t^) and macroscopic separations 
[Td <C Ts) differ from each other [5l|, [6y|- The different power law scahng of Td with respect to the separation q^ is 
illustrated by figure [TJ 



b? 0-1 ---^, - 




FIG. 1: Decoherence time scale Td for t^ ^ Ts at high temperatures in dependence of the separation \q — q'\ = 2go. For 
relatively small values of qo an thus Td > Tt ^ T^^, Td decays according to tV^ /{2'ykT)qQ (yellow dashed). If the decoherence 
time becomes very small Td <C tj,, the dependence ofTd on qo changes qualitatively to an inverse proportionality h / (y^ 'yF kT) q^ ^ . 
Here: F = lOwo (read), F = lOOtJo (blue). Further parameters: 7 = W~^mLOo, kT = bOhuJo- 



B. Long and short time behavior 



In this section we focus on the decoherence of macroscopic superpositions (large |a|) for the standard scenario 
Tb -^ Ts ^ Ty at high and low temperatures. The decoherence time scale Td thus has the same magnitude as the 
internal system dynamics Tg . For kT ^ huJa the functions (HO]) and (|151) become 



;,,^>,^,(i) = e-(4^'=^<'o/''^)* 



and 



Mrrf<r,(0 =e 



-{f'lkTql/h^)t 



(45) 



These approximated functions decay exponentially and differ just by a factor of two in the exponent. From figure 
UK one can see that the decoherence function ([55| follows the function iJ^TdKr^t) for short times t < Tg and for large 
times t ^ Ts oscillates around the approximation fJ-Td>TAt). The oscillations result from the rotation and breathing 
of the Wigner function in phase space (with frequency wo ~ ''"s^^) which is connected to a periodical change between 



superposition in coordinate and momentum space, 
become 



A^r,>r.W=e-(27-09.V^)* and ^lr,<rM 



For low temperatures kT <^ Hluo the functions (|40p and ([151) 



»ir^icn 2 



r 



7vh 



for i <C Tg , 
for t ^ Ts- 



(46) 



Figure [2b shows, that the decoherence function ([HS]) for short times t < Tg decay fast in accordance with the evolution 
of fJ-Td<Ts (Oi while the long-time behavior for i ^ r^ is again well approximated by fJ'Td>Ts (0- The long-time behavior 
of MTd<Ts(i) however can be characterized by a power law decay ^ t^^'^o , This result is in accordance with the 
behavior of a free quantum Brownian particle [67[. There is no characteristic decoherence time any more. In our case 
the decoherence function follows an algebraic decay for intermediate times and an exponential decay for large times. 




LJot 

{a.)kT = lOhujo, |aoP = 30. 




(b)fcT = 10--^?!t<;o, |aor = 100, 



FIG. 2: Time evolution of the decolierence function (I35|l (red line) in comparison to the approximation (I43p (blue dashed line) 
and (|40|) (yellow dashed line) for high temperatures (a) and low temperatures (b). The oscillations result from the rotation and 
breathing of the Wigner function in phase space (with frequency loo ~ t^^) which is related to a periodical change between 
superpositions in coordinate and momentum space. Further parameters: 7 = 10~^tJo, F — IOluq. 



C. Numerical Analysis of the dynamics 

1. Decoherence of a single-mode superposition 



Describing decoherence processes within the Born-Markov approximation is valid as long as the time scale of bath 
correlations Tb is the smallest time scale. Since the decoherence time scale r^ is indirect proportional to the bath 
temperature T and to the squared separation q^ , the decoherence process takes place on time scales that are comparable 
to the bath correlation time. In this case, non-Markovian influences become important. 

A gradual change from the Markovian to the non-Markovian regime can be studied by varying the bath correlation 
time in form of the inverse cutoff frequency F^^ and the decoherence time scale r^ by the parameter |ao| = 9o/2co. 
Such a change from the Markovian {t^ ^ ti,) to the non-Markovian regime {t^ ^ rt) is illustrated in figures [Stt-f by 
the evolution of fJ,i{t), fJ^Tdyrsit) ^nd IJ-TdKrA^)- Figures [3K-c show that a reduction of F slows down the decoherence 
process. For example, the decoherence time in fig. [3t for F = 2cl;o is five times larger than the decoherence time in 
fig. [3^ for F = lOOwo- More important are the qualitative changes in the evolution of the three decoherence measures 
fii{t), fJ-Tdyr^it) and t^Td<Tsit)- I^i figures[3K-d the norm /i/(t) follows ^-TdKrS'^), since the decoherence time is smaller 



than the characteristic system time scale Tg 
^ii{t) follows the evolution of /ird>r,(t). 



Wn 



A comparison to the case r^ 3> t^ is shown in figure [31i where 



So far, we have considered non-Markovian influences by approaching Td and r^, where the limiting cases for r^ 3> t^ 
and Td <C Ts have been distinguished. A further distinctive feature can be observed if additionally the condition 
Tb '3> Tg is fulfilled which means cjo ^ T- In this case the bath could be characterized as out-of-resonance [60]. The 
relaxation timescale t^ ^ jloq/T'^ also depends on the cut-off frequency but still remains by far the largest time scale. 





UJot 
(a)rd > n for T = IOOcjq, |ao| = 100. 



UJot 
(h)Td Si Ti for r = lOcJo, |ao| = 100. 
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UJot 
(c)rd Si r;, for T = 2<.Jo, |qo| = 100. 
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UJot 
(d)rd < Ti for r = 101.^0, |ao| = 500. 




(e)T6 <^Ts ^Ta for T = lOoJo, 
|aol = 10. 




UJot 



{i)n > Ts for r = O.OloJo, |aol = 30 

and 7 = O.OSloiq compared to 

7 = O.ltijo (dashed lines). 



FIG. 3: Time evolution of the decoherence function (I35|l (red lines) in comparison to the approximation (|43p (blue lines) and 
(|40[) fyeJJow lines) at high temperature kT = lOhujo and small coupling 7 — 10~^a;o in (a) to (e). Figures (a) to (d) illustrate 
the transition from the Markovian (rd 3> rt) to the non-Markovian regime (rd ^ Tb), related to the ratio of decoherence time 
Td and hath correlation time scale ti, « F^^. The decoherence function Hi{t) follows the approximated function (|43p since 
Td <C Tj, ~ ^o"^- A comparison with the case Td ^ Ts, where iJ,i{t) follows the function (|40ll is plotted in figure (e). Figure (/) 
illustrates the partial revivals of coherence in an out-of-resonance-bath. 



The time evolution of the decoherence function /i/(t) cannot be well approximated by ^J-Td<Ts or HTd>Ts alone. Figure 
|3f shows an example. For short times the norm fj.i{t) decays very fast, following the approximated function fJ.Td<T,- 
However, within the half of a system period ujgt ~ tt a partial revival of coherence takes place and the function fj.i{t) 
reaches a relative maximum that is given by the corresponding value of limit case fJ^T^yTs ■ Although the characteristic 
time scales of fJ,Td<T, and fJ,Td>Ts are quite different, only both limit cases taken together give a accurate picture of 
the decoherence process in this regime. The corresponding master equation (l27l) in this case is not of Lindblad-type 
with partially negative values of the coefficients 7|(i) and Ji{t). 



2. Disentanglement of a two-mode entangled coherent states 



In the following, we compare the findings on the decoherence process of a single-mode superposition state to the 
behavior of the concurrence of a two-mode entangled coherent state. Figures |4^-d show the time evolution of the 
concurrence Ci2(i) of an entangled coherent state |'1^+) = \a,d = 0,N = 2) in reservoirs with different specification 
of their non-Markovian character in form of the relation T/loq. Qualitatively the time evolution resembles strongly 
the results for a single-mode superposition cat state. Starting from a quasi Ohmic bath with (F ^ wq) in fig. UK) a 
reduction of the cut off frequency leads to deviations between the Markovian and non-Markovian results. For moderate 
non-Markovian influences the entanglement is preserved most efficiently (figHb). In the strongly non-Markovian out- 




1 1.5 2 

LOot 
(a)r = 3ujo, |ao| = 20. 





LOot 
(c)r = lO-^tjo, |ao| = 1500. 



(b)r = 0.5wo, |ao| = 20. 
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(d)r = 10~3(.j(,, laol =3500. 



FIG. 4: Time evolution of the concurreDce Ci2(t) (red line) for a two-mode entangled coherent state |"3?+) ((6| in a zero 
temperature environment with small ratios F/o^o and 7 = 10~^ci;o. Jn comparison the Markovian evolution is plotted for 
relaxation time 7"^^ (blue line) and with adjusted relaxation time r-y (yellow dashed line). In a moderate non-Markovian regime 
the entanglement is preserved for longer times (b) while it is lost faster in an out-of-resonance bath with F ^ ljo (c). The 
behavior of the concurrence is similar to that of the decoherence function. Revivals can occur even at zero temperature (d). 



of-resonance bath with T < luq this relation is only valid at short times while for longer times the non-Markovian 
concurrence even decays faster, as can be seen from fig. \^. The occurrence of coherence revivals is similar to the 
results found for the decoherence function of a single-mode superposition of coherent states. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the non-Markovian effects on decoherence and disentanglement processes of non- 
Gaussian continuous variable systems within the quantum Brownian motion model. We have compared the time 
evolution of the decoherence function of a single-mode cat state with the evolution of the concurrence of a two-mode 
entangled coherent state. For both cases we studied different decoherence and disentanglement scenarios depending 
on the relation between the characteristic time scales of system and environment. The entanglement dynamics of 
two-mode entangled coherent states is similar to the decoherence dynamics of single-mode coherent cat-states. We 
found exponential, Gaussian and algebraic decay patterns of the decoherence function in moderate non-Markovian 
regime and revivals of decoherence and concurrence in strongly non-Markovian out-of-resonance reservoirs. 
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